Breaking of the overall permutation symmetry in nonlinear optical susceptibilities of 
one-dimensional periodic dimerized Hiickel model 
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Based on one- dimensional single-electron infinite periodic models of trans-polyacetylene, we show 
analytically that the overall permutation symmetry of nonlinear optical susceptibilities is, although 
preserved in bound-state molecular systems, no longer generally held in periodic systems. The 
overall permutation symmetry breakdown provides a natural explanation of the widely observed 
large deviations of Kleinman symmetry in off-resonant regions of periodic systems. Finally, physical 
conditions to experimentally test the overall permutation symmetry breakdown are discussed. 

PACS: 78.66.Qn, 42.65.An, 72. 20. Dp, 78.20.Bh 



The nth-order optical susceptibility is generally defined as a rank n tensor x^d 1 a. 2 ---a n {— uji, 102, ■ • • ,Un), where 



Lo a — y u>i is the sum of incoming frequencies and \xa.\ ■ • • a n are the indices of spatial directions. The intrinsic 



permutation symmetry, as described in Butcher and Cotter's book [1], implies that the nth-order susceptibility is 
invariant under all n! permutations of pairs (ai,cJi), {0.%, a^),- • • , Intrinsic symmetry is a fundamental 

property of the nonlinear susceptibilities which arises from the principles of time invariance and causality, and applies 
universally to all physical systems. For the medium that is transparent and lossless for all relevant frequencies, i.e., 
far away from all transition frequencies, it is generally believed that the optical susceptibilities have a much more 
interesting property, namely, the overall permutation symmetry (or the full permutation symmetry in Boyd's book 
[2]), in which the susceptibilities are invariant when the permutation includes the additional pair Therefore, 
the nth-order susceptibility is invariant under all (n + 1)! permutations of the pairs (u, — u> a ), (ai, u>i), • • • , (a n , w n ) 
[3]. Furthermore, when the optical frequencies are much smaller than any of the transition frequencies, the dispersion 
of the medium at the relevant frequencies is negligible. It follows that the susceptibility is asymptotically invariant 
under all permutations of the subscripts fJ>, ct\, • ■ ■ , a n in the low frequency limit. This property is known as Kleinman 
symmetry [1,2,4]. 

However, as observed by Simpson and his coauthors [5], the overwhelming majority of recent optical experiments on 
organic materials and crystals showed large deviations from Kleinman symmetry, even in the low-frequency off-resonant 
regions. Although the deviations from Kleinman symmetry in molecular systems are fairly small. In this paper, based 
on the theoretical framework developed in our previous work [6-10], we prove that the overall permutation symmetry 
of nonlinear optical susceptibilities is broken in one-dimensional (ID) periodic systems. On the other hand, the overall 
permutation symmetry remains valid in bound-state molecular systems. Since the overall permutation symmetry is 
the basis of Kleinman symmetry, this provides a natural explanation why the deviations from Kleinman symmetry 
is much larger in periodic systems than in bound-state molecular systems. Indeed, despite the wide acceptance of 
the overall permutation symmetry in the nonlinear optics [1-3] , no direct measurement has tested the validity of the 
assertion. We will suggest physical conditions to experimentally test the overall permutation symmetry breakdown. 

The analytical derivations of the overall permutation symmetry of nonlinear susceptibilities are rigorous and correct 
in molecular systems [1-3] where the position operator r is well-defined in real space. However, for periodic systems, 
the usual definition of r is no longer valid over all space. Instead the "saw-liked" position operator must be introduced 
to maintain the periodic property of the system [11,12]. If periodic boundary conditions are applied to a physical 
system, the average electronic position could be anywhere for delocalized states [11]. This is clearly not the case for 
most molecular systems with only bound states. For periodic systems, the position operator r is conveniently defined 
in momentum space [13]: 
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r n k,n'k' = iVkCn.n' (k, k') +n nin -(k)5(k-k'), (1) 

where 

C„,„'(k,k')= / C^Wi'.k'W*. 
n„,„'(k) = - / <,k( r )VkMn',k(r)rfr, (2) 

V is the whole system volume, v is the unit cell volume, and V'n.kCr) = H n ,k(r)e !k ' r is the Bloch state with n and k 
the band index and crystal momentum, respectively. The two terms in Eq.(l) correspond to the intra- and inter-band 
transitions, respectively [13]. 

On the other hand, the p • A (current-current correlation) instead of the E • r (dipole-dipole correlation) gauge is 
often used, where p is treated as static [1,11]. Though quite successful in the linear case, the static current-current 
correlation is actually not equivalent to the dipole-dipole correlation (i.e., the E • r) gauge and gives the wrong results 
for periodic systems [8,9] in nonlinear optical studies. To restore the equivalence between these two gauges, one needs 
to incorporate the proper gauge phase factor in the current-current correlation [8] . The importance of the gauge phase 
factor in gauge equivalence is also verified in a recent work by Rzazcwski and Boyd [14]. For simplicity, we have used 
the E • r gauge in this paper. 

To demonstrate the symmetry breakdown, we focus on two ID single-electron periodic models - the Su-Shrieffer- 
Hccgcr(SSH) model [15] and its continuum sibling Takayama-Lin-Liu-Maki(TLM) model [16]. These two models are 
quite successful in explaining optical properties of trans-polyacetylene (see, for example, [17]). The advantages of ID 
single-electron periodic models are two folds. First, since £„ >n /(k, k') = 5 n ,n'5(k — k') can be applied in Eq.(l), these 
two models result in ID simple bands and are analytically solvable [18,19]. Second, single-electron models also satisfy 
physical conditions where the overall permutations symmetry is supposed to be strictly held in existing theories [1-3] , 
since the medium is loss-free at all the relevant optical frequencies under these models. When centro- or inversion 
symmetry is applied to the system, x^ vanishes and x^ becomes the first non-zero nonlinear susceptibility. Hence, 
we will concentrate our discussion on x ■ 

The Hamiltonian of the SSH model with the rigid-lattice approximation (i.e., the dimerized Hiickel model) is 
described as follows [15]: 



hssh — — 



*o + (-l)'| 



(cliA,s + clci + i,s), (3) 



where to is the transfer integral between the nearest-neighbor sites, A is the gap parameter and Cj s {Ci, s ) cre- 
ates (annihilates) a 7r electron at site I with spin s. Following the same procedure described in previous work [6-10], 
we consider the momentum space representation of the Hamiltonian given by Eq.(3). With the aid of the spinor 
description tpj, s (i)=(aJ, c s (£), o^" g (i)), where % s (t) and % s (t) are the excitations of electrons in the conduction band 
and the valence band with momentum k and spin s, we obtain the following formula: 

HssH(k, t) = Hg SH + H E r 

= E e{k)^l s {t)a^ ktS {t) - 1) ■ E e^\ (4) 

_ >L.<fc< » s 

2a — — 2a 1 

where 

d 

D = e (0( k )n,s^k,s+i-^Tpl s '>Pk,s), (5) 



e(k) — \J [2t cos(ka)] 2 + [Asin(ka)f , at (i = 1,2,3) are Pauli matrices, and f3{k) = —Atoa/e 2 (k). (3(k) is the 
coefficient related to the interband transition between the conduction and valence bands in a unit cell of length 2a, 
and the second term in Eq.(5) is related to the intraband transition [6-10]. 

According to field theory and the Wick theorem [20], the general four- wave- mixing(FWM) can be expressed as 
[6-10]: 
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x W ff (-n ;WllWa , WS ) = ^^ Yl J ^Trj(/3(fc)a 2+ *J,)G(fc, W ) 

9 1 

(P(k)<72 + ig^)G(k, ui - wi - w 2 - w 3 ) (6) 

where = u>i + iv 2 + <^3, L is the chain length, no is the number of chains per unit cross area, and V(u>i, uj 2 , uj 3 ) 
represents all six permutations over w\, lo 2 and oj 3 . The polymer chains are assumed to be oriented. Green's function 
G(k,u>) is defined as follows [6-8]: 

with = e(k)/h and e = + . 

The analytical expression for the general four- wave-mixing can be found in [10]: 

15 ^ f 1/s xdx 



V(z lt Z2,Z 3 ) ' 



y/(l-S 2 X 2 ){x 2 - 1) 



(2a; - -zi - z 3 ) 



x 8 (x - + z 2 )(£ - z 3 )(a; - z\ - z 2 - z 3 ) 

(2x + gi + z 3 ) ( 8 ) 

x 8 (x + zi)(x - z 2 )(x + z 3 )(x + z\ + z 2 + z 3 ) 

4(1 - 8 2 x 2 ){x 2 - 1) (3a; - 2^) (3a; - 2( Zl + z 2 + z 3 )) 

x 8 (x — z\ — z 2 ) (x — z\) 2 (x — z\— z 2 — z 3 ) 2 
4(1 - S 2 x 2 )(x 2 - 1) (3x + 2 Zl )(3x + 2(zi + z 2 + z 3 )) ) 
x 8 (x + Z! + z 2 ) (x + zi) 2 (x + z\ + z 2 + z 3 ) 2 y 



where x (3) ^ ^«!, Zj ^ (i = l, 2, 3), 5 = A and V{zu Z2iZ3) repreS ents all six permutations over z lt 

45 7r A" 2A 2^o 

22 and z 3 (Please note that the summation over all six permutations was missing in the original Eq. (3.3) in Ref. [10] 
due to a typo.). 

To discuss the breaking of the overall permutation symmetry, we may simply consider the difference 

<5x (3) =Xs 3 SH(- rj ; a, i' w 2,w 3 )-xi5H(wi;-ri,w 2 ,w 3 ). (9) 

Due to the trivial permutation of the spatial indices, the spatial indices fj,a\ ■ ■ ■ a n can be dropped in the expression 
of x (n) in the ID case. Thus, the Kleinman symmetry is exactly the same as the overall permutation symmetry for 
the ID case, though it is generally different from the overall permutation symmetry for higher dimensional cases. 
Obviously, <5x' 3 ^ in Eq.(9) should be identically zero if the overall permutation symmetry were preserved. However, it 

is easy to see from (8) that Xssh(~ ^> Wl ' W2 > w 3) ^ s symmetric in u)\, u 2 , and ui 3 , while Xssir^ 1 ' — ^1^2,^3) is only 
symmetric in uj 2 and lu 3 . Hence, can not be zero everywhere. Indeed, replacing oj\ by — (o>i + uj 2 + uj 3 ) in (8), 
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substituting the resulting expression into (9), and then simplifying the final expression, we obtain 

<$X (3) = XsSffM^l + ^2 +UJ 3 )]UJi,LU 2 ,LU 3 ) - Xsiff^i; -(<*>! +^2 + W 3 ),W2,W 3 ) 



(3) J15_ ^ f 1 ^ V(l - 8 2 x 2 )(x 2 - l)da; 



Xo 256 



(3x + 2z 2 )(3x + 2(zi + z 2 + z 3 )) / 1 1 



(x + z 2 ) 2 (x + z\ + z 2 + z 3 ) 2 \X + Zl + Z 2 x + z 2 + z 3/ 
(3x - 2z 2 )(3x - 2(zi + z 2 + z 3 )) / 1 1 \ 



(a; - z 2 ) 2 (a; - zi - z 2 - z 3 ) 2 \a;-zi-z 2 a; - z 2 - z 3 
(3x + 2zi)(3x-2z 2 ) ^ 1 1 



(x + zi) 2 (x - z 2 ) 2 \x + zi + z 3 x-z 2 -z 3 
(3x-2zi)(3x + 2z 2 ) /ll 



(x - zi) 2 (x + z 2 ) 2 \x — Zl — z 3 X + z 2 + z 3 

where the summation is over the two permutations of z 2 and z 3 . Eq. (10) provides a general expression for the 
deviation of the overall permutation symmetry under the SSH or ID dimerized Hiickel model. 

Remark .1. A closer examination shows that the first two terms in (8) come from the interband transition, and the 
last two terms in (8) come from the intraband transition. The deviation of the overall permutation symmetry (Eq. 
(10)) is due to the intraband transition alone. 

To demonstrate the large deviation of the overall permutation symmetry quantitatively, we may consider two special 
cases under the TLM model (the continuum limit of the SSH model): Xtlm (~ 3c; oj, oj, oj) and Xtlm^'i w ' w ' — 3c). 
Our analysis shows that 



(11) 



and 



( 3) , 5 80 2 98580 4 n . , 



X ^ 3 | M (c; c, c, -3c) = x < 3 ) -1^1^(40^ - 61)/(z) + ^(4z 2 - l)/(2z) 

-^(1944z 2 - 241)/(3z) + ||(27z 4 - 30z 2 + 805) J < 12 > 

(3 ), 5 80 , 28500 , fixx , 

= ^<28 + 33* + TooT* + ( ^ 0) - 

where z = fic/2A, and the function / is defined by the formula 

arcsin(z) 

/(*) = < 



(z 2 <l) 

cosh (z) VK , v 7 



2; 



\/I 2— I 2zV~z 



For polyacetylene, using typical parameters of to = 2.5eV, A = 0.9eV, no — 3.2 x 10 14 cm~ 2 and a = \.22A results in 
~ 1-0 x 10~ 10 esu [6-10]. In Fig.l we have plotted Xtlm(~ 3w; c,c,c), Xtlm^' u,uj, — 3c) and the difference 
between the above two x^ m the off-resonant region. The graph shows that there is about a 40% difference for 
0=1/6 (about 0.3eV or 4.14/xm) between these two quantities. Subtracting the two asymptotic expressions in Eq. 

(11) and (12) reveals that the difference between Xtlm(~ ^ lu '- w ' w ' w ) anc ^ Xtlm( u > w ' w ' — 3c) in the off-resonant 
region satifies the following relationship: 

^)M« e4n °y v . (M) 
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We have also computed nonlinear susceptibilities without the Vfc or dk terms (corresponding to the intraband 
transitions) in Eq.(5) [10]. The results preserve the overall permutation symmetry. Excluding the gradient terms is a 
crude approximation for molecular systems where the polarization current is not present [21]. Therefore in periodic 
systems, it is the gradient term that breaks the overall permutation symmetry which remains valid in molecular 
systems [1-3]. Obviously, this is closely related to the fact that the position operator for periodic systems is entirely 
different from that for molecular systems. 

Remark .2. Recently, deviations from Kleinman symmetry in the low- frequency off-resonant regions have been ob- 
served in many nonlinear optical experiments [22-31]. Upon a careful examination of these experiments, we observed 
that: (i). The deviation of the Kleinman symmetry increases with decreasing band gap and is proportional to oJ 2 for 
crystals [22]; (ii). the deviation of the Kleinman symmetry in dclocalized states such as aromatic molecules [28,29] 
and some polymers [25-27] or crystals [22-24] is usually much larger than that in localized states such as molecular 
systems such as O2, N 2 , etc [30,31] (20^50% versus <8%). Eq.(14) captures the dependence of uj 2 and the bad gap 
observed in (i). Furthermore, since the overall permutation symmetry is the basis of the Kleinman symmetry, our 
computation also provides a natural explanation about (ii), as we have shown that the overall permutation symmetry 
is broken in periodic systems. Finally, the vanishing x^ under the SSH or TLM model shows that some symmetries 
such as centro-symmetry can suppress the deviation from Kleinman symmetry even for periodic systems. This may 
explain why Kleinman symmetry is still preserved in some x^ experiments of crystals [32]. 

Kleinman symmetry (or Kleinman conjecture) is considered valid only when u> = [1,2]. Previous theoretical 
calculations [33-35] have explicitly shown the Kleinman symmetry breakdown when u> 7^ for higher dimensional cases. 
Nevertheless, it is still generally treated as a valid approximate symmetry for all physical systems in low-frequency 
off-resonant regions. However, many experiments show that the range of Kleinman symmetry is substantially more 
restrictive than what is widely assumed in theory, especially for periodic systems [5]. To explain the large Kleinman 
symmetry deviation in periodic systems [22-31], various models have been presented. For example, Levine's model 
that predicts the second order polarizability tensor d F /d A cx J 1 where superscripts F and A mean forbidden and 
allowed respectively [36] , dipole contributions from two perpendicular directions [28] , harmonic frequency 2w strongly 
resonant with the Q band [29] , mutual exclusion properties between Kramers-Kronig dispersion relations and Kleinman 
symmetry [5], etc. Without so many input parameters, the overall permutation symmetry breakdown in periodic 
systems provides a much more straightforward and general explanation of the above experiments. 

In conclusion, the overall permutation symmetry for nonlinear susceptibilities is, albeit preserved in bound-state 
molecular systems, no longer generally held in periodic systems with delocalized states. Therefore, it leads to large 
deviations of Kleinman symmetry. The deviation Sx^ is proportional to (e^not^^h 2 ^ 2 )/ A 8 in the off-resonant 
regions. Theoretically, non-interacting centro-symmetric ID periodic structures such as single crystals of conjugated 
molecules are ideal materials for performing off-resonant x^ experiments to directly test the overall permutation 
symmetry breakdown. Practically, electron-electron correlation and other interactions are very important in ID 
conjugated systems [17]. These interactions will certainly change the magnitude of the deviation. However, as long 
as those interactions are not strong enough to destroy the band structures where the intraband transitions remain 
valid, the overall permutation symmetry breakdown can still be observed. 
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